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Using nuclear quadrupole resonance, the phase diagram of 1111 7?FeAsOi_a:Fa; (i?=La, Ce, Sm) iron pnic¬ 
tides is constructed as a function of the local charge distribution in the paramagnetic state, which features 
low-doping-like (LD-like) and high-doping-like (HD-like) regions. Compounds based on magnetic rare earths 
(Ce, Sm) display a unified behavior, and comparison with La-based compounds reveals the detrimental role 
of static iron Sd magnetism on superconductivity, as well as a qualitatively different evolution of the latter at 
high doping. It is found that the LD-like regions fully account for the orthorhombicity of the system, and are 
thus the origin of any static iron magnetism. Orthorhombicity and static magnetism are not hindered by su¬ 
perconductivity but limited by dilution effects, in agreement with 2D (respectively 3D) nearest-neighbor square 
lattice site percolation when the rare earth is nonmagnetic (respectively magnetic). The LD-like regions are not 
intrinsically supportive of superconductivity, on the contrary of the HD-like regions, as evidenced by the well- 
defined Uemura relation between the superconducting transition temperature and the superfiuid density when 
accounting for the proximity effect. This leads us to propose a complete description of the interplay of ground 
states in 1111 pnictides, where nanoscopic regions compete to establish the ground state through suppression of 
superconductivity by static magnetism, and extension of superconductivity by proximity effect. 

PACS numbers: 74.70.Xa, 75.25.Dk, 76.60.-k 


I. INTRODUCTION 

One of the main motivations to study iron-based high- 
temperature superconductors (IBS) is the frequent vicinity 
of unconventional superconductivity and static magnetism in 
their phase diagram. Also observed in copper-based super¬ 
conductors, it raises the questions of ground state interplay 
and of the role of spin fluctuations in the Cooper pairing in¬ 
teraction. By applying pressure, changing the doping or alter¬ 
ing the electronic structure via isovalent substitutions [1, 2], 
the ground state of most IBS can be modified from static 
magnetism towards superconductivity, with the possibility of 
ground state coexistence on a microscopic scale in several 
cases [3-10]. 

A key aspect of iron-based superconductors is the itiner¬ 
ant, multi-band character of their electronic structure. In 
compounds such as the 1111 and 122 families, the low- 
temperature magnetic order of the parent compounds (T^v « 
140 K for the 1111 family) is argued to be a spin-density 
wave caused by the nesting properties of the Fermi surface 
[11, 12]. As this transition is associated to a tetragonal-to- 
orthorhombic structural transition at the same or a slightly- 
higher temperature, both Ising-nematic spin fluctuations and 
orbital fluctuations have been invoked to explain the observed 
anisotropic electronic response of the iron plane above the 
magnetic transition [13]. Evolutions of the Fermi surface with 
doping would affect the competition of ground states, and rel¬ 
ative changes in interband and intraband scattering would play 
a role in determining whether superconductivity is associated 
to spin or orbital fluctuations, yielding respectively a s± or 
S++ symmetry of the order parameter [14—17]. 

The ground state competition is, however, complicated by 
spatial electronic inhomogeneities which derive either from 


the local effect of in-plane chemical substitutions [18-20], or 
from intrinsic iron plane physics. Using Nuclear Quadrupole 
Resonance (NQR), a local probe of the charge environment, 
we have previously shown that the fluorine-doped 1111 fam¬ 
ily (i?FeAsOi_a;Fa; with R a rare earth) features two types 
of local electronic environments intrinsic to the FeAs layers 
when moderately doped [21]. Although these environments 
are already defined in the paramagnetic state, their simulta¬ 
neous presence in purely magnetic or purely superconducting 
samples appears to rule out a trivial explanation of any ground 
state coexistence in 1111 compounds. Further complicating 
the analysis of ground state interplay, such coexistence seems 
dependent on the choice of rare earth (La, Ce, Sm...) [3, 22]. 
Obscuring the whole phase diagram is the difficulty in assess¬ 
ing the true doping, which may differ dramatically from the 
nominal doping depending on the rare earth and the synthesis 
route [23-25]. 

In this article, we present how the nanoscale charge envi¬ 
ronments observed by NQR in the paramagnetic state yield 
a highly-accurate determination of the effective doping, al¬ 
lowing us to compare unambiguously the phase diagrams of 
i?FeAsOi_a;Fa; (i?=La, Ce, Sm) and to explain the evolu¬ 
tion of ground state properties with doping. Sec. Ill describes 
how to obtain such phase diagrams, which are found undis- 
tinguishable for the two magnetic rare earths (Ce, Sm). The 
different behavior for lanthanum shows that static magnetism 
from the iron has a detrimental effect on superconductivity, 
while at high doping La-based samples display a specific, non- 
monotonous behavior of the superconducting transition tem¬ 
perature. In Sec. IV, it is found that the low-doping-like re¬ 
gions seen by NQR fully account for the orthorhombicity and 
iron static magnetism of the system. Nearest-neighbor square 
lattice percolation is at play, with dimensionality being depen- 
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dent on whether the rare earth itself is magnetic. In Sec. V, 
it is shown that the low-doping-like regions can only host su¬ 
perconductivity by proximity, with intrinsic superconductivity 
originating from the high-doping-like regions seen by NQR. 
Finally, we propose in Sec. VI a full description of the in¬ 
terplay of ground states in 1111 pnictides, where supercon¬ 
ductivity extends spatially whenever static magnetism is weak 
enough in the nearby microscopic regions. 


II. EXPERIMENTAL DETAILS 


Two different routes were followed for the synthesis of the 
polycrystalline i?FeAsOi_a;Fa; (i?=La, Ce, Sm) samples. For 
route 1, FeAs was prepared by a solid-state reaction prior to 
the synthesis of the corresponding 1111 compounds. Subse¬ 
quently, the resulting FeAs was mixed with metallic R, R 2 O 3 , 
and RF 3 in a stoichiometric ratio. This mixture was homoge¬ 
nized by grinding in a mortar. For route 2, we prepared RAs 
as first step by reacting R and As lumps in a stoichiometric 
ratio via a vapor transport reaction. The second step of route 
2 used the resulting RAs, Fe, Fe 203 , and FeFa as starting 
materials in a stoichiometric ratio. Here, the starting mate¬ 
rials were homogenized by grinding in a ball mill. In either 
case, the resulting powders were pressed into pellets under 
Ar atmosphere, and subsequently annealed in an evacuated 
quartz tube in a two-step synthesis first at 940° C for 12 h 
and then at 1150°C for up to 60 h. All samples were charac¬ 
terized by powder X-ray diffraction (XRD) and by scanning 
electron microscopy (SEM) with semiquantitative elemental 
analysis using the wavelength dispersive X-ray (WDX) mode. 
Their structural properties and their structural, magnetic, and 
superconducting phase transitions were determined using a 
wide array of techniques: standard XRD, synchrotron XRD at 
low temperature, magnetic susceptibility, electrical resistivity, 
muon spin relaxation and rotation (pSR), Mossbauer spec¬ 
troscopy, and nuclear magnetic resonance (NMR) [22, 26- 
30]. The determination of the doping will be discussed in 
the next Section, with all phase transitions being included in 
Fig. 3. 

In order to probe the distribution of charge environments 
at the nanometer scale, ^®As nuclear quadrupole resonance 
(NQR) was used. NQR probes nuclei with spin I> 1/2, which 
possess a finite electric quadrupole moment eQ with e the el¬ 
ementary charge and Q the nuclear quadrupole moment. In 
the absence of applied or internal magnetic field, any finite 
electric field gradient (EFG) at the atomic site will lift at least 
partially the degeneracy of the nuclear electric quadrupole en¬ 
ergy levels, according to the following Hamiltonian [31]: 


eQV.z 
4/(2/ - 1) 



/(/ + l) + |(/2 +/2) 


( 1 ) 


whose coordinate axes are defined by the principal axes of 
the electric field gradient (EEG) tensor, itself characterized by 
its largest principal value Vzz and its asymmetry parameter 
r]={V.. - Vyy)/Vzz (|E,,| > \Vyy\ > |V,,| and 0 < 77 < 1 ). 
Since ^®As has a nuclear spin /=3/2, the Hamiltonian yields 


a single resonance frequency: 




eQVzz 

2 h 
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where h is Planck’s constant. The electric quadrupole fre¬ 
quency i/Q depends on the symmetry and amplitude of the 
EEG, which itself depends on the atomic/nanoscale charge en¬ 
vironment. Eor every local charge environment as seen from 
the arsenic sites, one resonance line will be observed. 

Using a standard pulsed NQR spectrometer, radiofrequency 
irradiations over a typical ^'r/=8-14 MHz frequency range 
were performed to obtain the histogram of resonance frequen¬ 
cies for all arsenic sites in the sample. Powder samples were 
crushed to ensure that the crystallites are small enough to ob¬ 
tain a good penetration of the radiofrequency field. Note that 
operating on powders does not affect the spectra, since R de¬ 
pends on the principal axes of the local EEG, i.e., the reso¬ 
nance frequency for a given charge environment is indepen¬ 
dent of the orientation of individual crystallites. Spin echo se¬ 
quences (|-r-7r) were used, with typical t=20-30 /is. The t 
values and the pulse sequence repetition rates were checked to 
be small enough that spin-spin or spin-lattice relaxation con¬ 
trast do not distort the relative line intensities. Point-by-point 
spectra were obtained by integrating the full echo and apply¬ 
ing a intensity correction. 


III. RECONSTRUCTED PHASE DIAGRAM 


In most of the published literature, nominal dopings are 
used to build the phase diagrams of iron-based superconduc¬ 
tors. While this yields qualitatively correct results for a given 
sample series, it is known that the real doping may be signifi¬ 
cantly lower, especially for high nominal dopings [23-25]. It 
is thus difficult to extract quantitative information (e.g.: dop¬ 
ing thresholds), and to compare different compounds in terms 
of ground-state competition and the doping profile of the su¬ 
perconductivity dome. Eor E-doped 1111 compounds, real 
content determination using WDX spectroscopy is hardly pos¬ 
sible for cerium-based compounds due to the superposition of 
the relevant Ce and E lines, and may generally be affected 
by the presence of poorly-crystallized impurities which XRD 
cannot easily account for. Using fluorine NMR, the doping 
could be established within an absolute error of 2% [6, 32], 
which is satisfactory for high dopings but may be imprecise 
to discuss the boundary between the two ground states at low 
doping. Here, we propose to take advantage of the doping de¬ 
pendence of the charge environments previously seen by NQR 
for La- and Sm-based samples [21]. 


A. Principle 

Eig. 1 shows the NQR spectra for /?-llll (R = La, Ce, 
Sm) in the paramagnetic state over a broad doping range (it¬ 
self determined by NQR as described below, see also Tab. I 
for the fluorine content). Whereas undoped and highly-doped 
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Figure 1: (Color online) NQR spectra of the studied samples in the paramagnetic state. The doping is indicated as xnqr/xfswt (FSWT: 
full spectral weight transfer, see text). The measurement temperature is coded by the color (see also the scale in Fig. 2): Yellow for room 
temperature, magenta for 150 K to 160 K, indigo for 100 K (Ce, xnqr/xfswt=0.63), blue for 40 K to 55 K. Full lines are Gaussian fits. 
Some low-temperature La/Sm spectra are from a previous study [21]. 


samples display a single peak corresponding to a single charge 
environment, the samples at intermediate doping display two 
peaks (two charge environments) which were shown to cor¬ 
respond to a nanoscale electronic inhomogeneity intrinsic to 
the FeAs layers [21]. Defining Fq {vq) and (w^) as the 
low-frequency (high-frequency) peak’s position and spectral 
weight {w^+w^=1), and x as the WDX-established doping, 
these two charge environments were found to have the follow¬ 
ing properties: (i) the transfer of spectral weight is linear with 
X and is completed {w^ = 0) at a:«0.1 (ii) Vq and Vq vary lin¬ 
early with X, extrapolating respectively to the frequencies of 
the x=0 and a;=0.1 charge environments (iii) the linear x de¬ 


pendence of Fq appeared to extend beyond a;=0.1, which was 
independently confirmed [33]. The low- and high-frequency 
peaks were thus associated to low-doping-like (LD-like) and 
high-doping-like (HD-like) regions. 

Taking advantage of the fact that NQR is immune to spu¬ 
rious phases and allows to measure with good precision the 
peak positions and the spectral weights in practically all sam¬ 
ples, we reverse the above observations to define xnqr, a 
NQR-defined doping. For samples with doping low-enough 
that two charge environments are observed, following point 
(i) we set: 

Xnqr = w^xrswt. 


( 3 ) 
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Figure 2: (Color online) Quadrupolar frequency versus xnqr/xfswt ■ The L/H 1 /H 2 /H 3 labels refer to the observed frequency branches. 
Literature data is taken from Refs. 33-38 (for one sample, the unknown measurement temperature is color-coded as white). For 
Xnqr/xfswt<^, xnqr is determined using the NQR spectral weights (see text for comments about some samples from Ref. 33). Linear 
fits are indicated as full lines with attached labels. For each sample, the error on xnqr is derived from fits of the spectra at one or more 
temperatures. For xnqr/xfswt>^, xnqr is determined from the linear extrapolation of the low xnqr behavior, indicated as a dashed 
line. For each sample, the error on xnqr is set so as to account for the uncertainty on rq, whether experimental on a single point or due to 
remaining temperature dependence. 


with xpswT the value for which full spectral weight transfer 
(FSWT) has occurred, i.e., the electronic homogeneity is re¬ 
stored. For samples with doping large enough that a single 
charge environment is observed {w^=\), following point (iii) 
we set Xnqr such that Rq is linear versus this variable. 

While it seems from point (i) that x fswt corresponds to 
about 10 % fluorine content, we will show in the following 
that this value cannot be estimated very precisely and that it 
is more reliable to express xnqr in units of xrswt (he., to 
determine the ratio Xnqr/xrswt) rather than as a fluorine- 
content equivalent doping. 


B. Data extraction 

From all measured spectra and using our rules for deter¬ 
mining Xnqr/Xfswt, we obtain Fig. 2, where the verti¬ 
cal axes are offset so that there is overlap for all three rare 
earths at Xnqr=Q and room temperature. Such overlap is pre¬ 


served for the whole Vq low-frequency branch (“L”) at room 
temperature, with a well-defined linear dependence on x nqr 
(yellow fits) as expected from the above points (i) and (ii). 
This linear behavior holds against changes in temperature, al¬ 
beit with different slopes for the different rare earths (pur¬ 
ple fits for r«160 K). The influence of temperature will be 
considered when discussing orthorhombicity and static mag¬ 
netism in Sec. IV. For the Vq high-frequency branch with 
xnqr<xfswt, the situation is complicated by the pres¬ 
ence of a superstructure in the high-frequency peak of the 
Sm samples, which appears to be partially echoed in the 
xnqr/xfswt=0-55 La sample and xnqr/xfswt=0-63 
Ce sample, thus a total of three high-frequency branches 
(“H1/H2/H3”). Nonetheless, similar linear behaviors account 
well for the data (see yellow and purple fits). For the sam¬ 
ples with xnqr>xfswt, he., only La- or Ce- based sam¬ 
ples, we apply the above procedure of extrapolating the lin¬ 
ear behavior at xnqr<xfswt- The Hi frequency branch is 
used, since the H 2 and H 3 branches appear to play little to no 
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R Source 

X 

xnqr/xfswt 

Xfswt 

La this study 

0 

0.00 



0.035 

0.31(4) 

0.11 


0.04 

0.37(4) 

0.11 


0.045 

0.41(2) 

0.11 


0.05 

0.49(4) 

0.10 


0.075 

0.55(3) 

0.14 


0.1 

1.19(6) 

0.08 


0.15 

1.90(2) 

0.08 

Ce this study 

0 

0.00 



0.05 

0.43(2) 

0.12 


0.15 

0.50(2) 

0.30 


0.1 

0.63(2) 

0.16 


0.1 

0.95 

0.11 


0.2 

1.94 

0.10 


0.25 

2.12(7) 

0.12 

Sm this study 

0 

0.00 



0.04 

0.34(2) 

0.12 


0.06 

0.59(2) 

0.10 


0.08 

0.65(4) 

0.12 


0.1 

0.90(4) 

0.11 

La [34] - Kitagawa 

0.14 

1.98 

0.07 

[33] - Oka 

0.03 

0.40(4) 

0.07 


0.04 

0.48(4) 

0.08 


0.06 

0.90(5) 

0.07 


0.08 

1.09 

0.07 


0.1 

1.67 

0.06 


0.15 

2.64(8) 

0.06 

[35] - Tatsumi 

0.14 

1.51 

0.09 

[36] - Mukuda 

0.22 

2.03(22) 

0.11 

[37] - Kobayashi 

0.11 

0.56(5) 

0.20 


0.15 

1.33 

0.11 

[38] - Nitta 

0.14 

1.80(27) 

0.08 

Table I: Doping information 

on the samples from this study, as well 


as on samples from previous LaOi_a;Fa;FeAs NQR/NMR studies 
whose spectral properties have been re-analyzed following the same 
procedure as for our samples. R indicates the rare earth. For samples 
from the literature, the name of the first author is given next to the 
reference number, x indicates the fluorine doping. For our La-based 
samples it is the nominal value, which is on the order of the WDX 
value, whereas for our Sm-based samples it is a rescaled value de¬ 
rived from WDX measurements [26, 27]. For our Ce samples it is the 
nominal value, without the possibility of comparison to WDX values 
due to the superposition of the relevant Ce and F lines. For samples 
from the literature it is usually the nominal value, xmqr/xfswt is 
the NQR-derived doping as defined in Sec. Ill A. xfswt is the per- 
sample extrapolated fluorine content for which full spectral weight 
transfer would occur in the NQR spectrum (see Sec. IIIC). 


role for La- and Ce-based samples. Note that this may not 
hold true for Sm-based samples (no sufficiently highly-doped 
samples were available to test for this). Since at high doping 
the temperature dependence of the spectra vanishes (La) or is 
much reduced (Ce), we use a single extrapolation intermediate 
between the low-temperature (La, Ce) and room-temperature 
(La) Hi fits (no reliable room-temperature Hi fit is available 
for Ce, but the corresponding data is roughly similar to the La 
data). This interpolation is shown as a dashed line on Fig. 2. 

The procedure to determine xmqr/xpswt was also ap¬ 
plied to LaOi_a;Fa;FeAs NQR/NMR data from the literature 
(Refs. 33-38). All included samples are listed in Tab. I to¬ 
gether with their Xnqr/xfswt value. In the following, they 
are designated as “source-x”, where source is the first author 
of the original paper and x the reported fluorine content. In 
most cases, the availability of the NQR spectra allowed us to 
re-analyse them in the same way as for our own samples. A 
few samples for which data was too sparse or fit anomalies 
were present are excluded. For the samples with two-peaked 
spectra (Oka-0.03, Oka-0.04, and Kobayashi-0.11), the spec¬ 
tral weights obtained from the fits yield xnqr/xrswt values 
respectively equal to 0.40(4), 0.57(2), and 0.56(5). However, 
the reliability of these values may be affected by the exper¬ 
imental conditions under which the spectra were measured, 
i.e., the radiofrequency pulse sequence repetition time and the 
delay r between the pulses (see Sec. II). If the former is too 
short or the latter too long, the relative spectral weights will be 
modified due to differences in Ti (spin-lattice relaxation time) 
or T 2 (spin-spin relaxation time) between the two peaks [21]. 
This may be the case for the Oka-0.06 sample, whose single 
NQR peak’s frequency implies xnqr/xrswt—0-90(5), i.e., 
a value for which two peaks should be observed. Here, the 
expected low-frequency peak may be entirely missed due to 
its small weight (10(5)%) being further reduced by so-called 
T 2 relaxation contrast. Considering the resonance frequencies 
obtained from the fits, the above Xnqr/xrswt values for 
Oka-0.03 and Kobayashi-0.11 appear reasonably compatible 
with the L and Hi frequency branches defined by our sam¬ 
ples (see the corresponding diamond symbols on Fig. 2). For 
Oka-0.04, the agreement is somewhat poorer and suggests that 
the actual Xnqr/xrswt value is lower than the fit-derived 
value of 0.57(2). This is supported by the comparison of spin- 
lattice relaxation rate measurements in Ref. 33 and Ref. 21, 
showing that Oka-0.04 is indeed rather lower doped than 
our 5% fluorine sample for which xmqr/xfswt=0-A9(A). 
Based on agreement with the L and Hi frequency branches, 
xnqr/xfswt is then estimated to be 0.48(4), as used in 
Fig. 2. 


C. Estimation of a;FSWT 

All obtained xmqr/ xpswT values are reported in 
Tab. I. To get an estimate of xfswt, one can compute 
x(xnqr/xfswt)~^ for every sample; by setting xnqr=x, 
Xfswt will be the fluorine content for which full NQR spec¬ 
tral weight transfer would occur in a given sample, with an 
inaccuracy corresponding to that present in x. As can be seen 
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in Tab. I most computed values of xpswT are in the range 
0.06-0.12, while a few values are slightly (0.14, 0.16) or much 
(0.20, 0.30) higher. The most likely candidate for this variabil¬ 
ity is fluorine that did not enter the matrix 1111 phase (push¬ 
ing xpswT up for a given sample), although other sources 
cannot be ruled out (e.g.: small amounts of oxygen or arsenic 
vacancies which might drive down xfswt)- This variability 
plays however no role in our study, where all deductions will 
be made on the basis of the Xnqr/xfswt ratio, which does 
not depend on any chemical characterization. 


D. Phase diagram 

The phase diagram for all three rare earths is rebuilt in a 
systematic way as a function of Xnqr/xfswt, as shown 
on Fig. 3. Tc data for La-based samples have been scaled 
so that the maximal transition temperature TJ"““(La)=26 K 
coincides with T™“^(Sm)=55 K [39]. For cerium, the scal¬ 
ing of Tc data is done so that there is good overlap with the 
data for Sm-based samples. The characterization of our sam¬ 
ples is taken almost exclusively from references [22, 26-30], 
in particular the work of Maeter et al. where hrst estimates of 
xnqr were used for Ce-based samples. For samples from the 
literature, the transition temperatures are taken from the orig¬ 
inal papers. For a given technique, the data was harmonized 
whenever possible by extracting the transition temperature in 
a systematic way (e.g.; midpoint of the resistivity drop for the 
superconducting transition). In any case, the use of different 
techniques for a given phase transition necessarily yields some 
additional point dispersion, which is taken into account by the 
visual guides of Fig. 3. The transition lines are well-defined, 
which shows the adequacy of our NQR approach to the dop¬ 
ing. This is especially true for La-based compounds, with 
twelve samples from other NQR groups with different sam¬ 
ples sources being successfully combined to our eight sam¬ 
ples. 

The tetragonal-to-oithorhombic structural transition and 
the transition to long-range-ordered (LRO) magnetism can 
be described by single Ts and T^i^lro lines for cerium and 
samarium. In-between these two transitions, a single Tm,sro 
line describes the transition to short-range-ordered (SRO) 
magnetism, which extends slightly into the superconductiv¬ 
ity region. A good scaling of the Ce and Sm Tc data can be 
achieved, with a scaling ratio « 0.73 which is similar to the 
ratio of the maximal transition temperatures (43/55 « 0.78) 
[29, 39]. While no data is available for highly-doped Sm- 
based samples, this is evidence that the superconducting re¬ 
gions for Ce and Sm are likely very similar, with maximal Tc 
being reached only at dopings far above the range where static 
magnetism and superconductivity meet. Finally, static mag¬ 
netism of the rare earths also extends far beyond this range. 
Thus, the phase diagrams for these two magnetic rare earths 
show no sign of differences after rescaling Tc- 

In the case of lanthanum, several differences occur. 
Tn.lro would appear to decrease faster: while few sam¬ 
ples show this transition, there are several low-doped super¬ 
conducting samples which do not show it, therefore con¬ 


straining it to below xnqr/xfswt^O-^- SRO magnetism 
is completely absent from superconducting samples, although 
it cannot be ruled out in magnetic samples as a precursor 
of the LRO magnetic transition. Finally, the superconduc¬ 
tivity region starts at lower doping and grows only until 
Xnqr!xfswt~^, after which it decreases with a sign of 
recovery only at very high doping. A candidate for the dif¬ 
ference in static magnetism compared to the unified behav¬ 
ior of Ce- and Sm-based samples is the lack of magnetism of 
lanthanum, whose formal 3 h- charge implies closed electronic 
shells. This is verified experimentally by spin-lattice relax¬ 
ation rate measurements: the spin fluctuations in La-based 
samples are much smaller (as well as more doping-dependent) 
than in samples containing a magnetic rare earth (Ce, Sm, 
Nd), and are thus ascribed to iron plane magnetism for the 
former and to 4/ rare earth magnetism for the latter [21, 40- 
43]. Static magnetism in the iron planes may thus be promoted 
at higher doping by magnetic rare earths (to be discussed in 
Sec. IV), and the concomitant delaying of superconductivity 
is a strong indication that static magnetism is detrimental to 
superconductivity. Regarding superconductivity at high fluo¬ 
rine doping in La-based samples, while it was already known 
that Tc tends to decrease [27], our phase diagram shows to our 
knowledge the first indication of eventual Tc recovery. This 
is compatible with observations on hydrogen-doped samples 
[44]: a second superconductivity dome occurs for La at high 
doping, whereas for Ce and Sm there is a single dome extend¬ 
ing to more than 40% doping. This difference is then argued 
to be related to the doping-dependence of the iron id orbital 
degeneracy. 

Note that the “ribbons” of Fig. 3 should be seen as an 
interval of confidence for each transition line. Especially, 
the T—^0 superposition of the magnetic and superconducting 
phase transitions about xnqr/xfswt^OA (La) and «0.5 
(Ce, Sm) should not be understood as a suggestion of a quan¬ 
tum critical point, whose relevance is anyhow unclear in an 
inhomogeneous context. There may be a crossing of the two 
transitions lines, which must however be very limited in dop¬ 
ing range considering that no sample is known to feature both 
long-range magnetic order and superconductivity. This will 
be discussed in terms of ground state competition in Sec. VI. 

Considering the electronic separation into regions which 
seem related either to the magnetically-ordered, undoped 
compound or to the superconducting-only Xmqr/Xfswt=^ 
compound (see Fig. 2), it is a distinct possibility that these re¬ 
gions are respectively more prone to develop static magnetism 
or superconductivity. Taking into account our observations 
about the phase diagram and the fact that most samples ex¬ 
hibit only one of the two ground states, we examine this issue 
in the next two sections. 


IV. ORTHORHOMBICITY AND STATIC MAGNETISM 
A. Spatial origin 

As shown on Figs. 1 and 2, the spectra in the param¬ 
agnetic state display a sizeable temperature dependence for 
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Figure 3: (Color online) Phase diagram of 7?FeAsOi_a:Fa; (i?=La, Ce, Sm) as a function of xnqr/xfswt ■ Tc data for La- and Ce-based 
samples have been scaled (see text). For some samples and some transitions, more than one point may be displayed due to characterization 
using multiple techniques. The samples from the literature are referenced in Tab. I, except the undoped Ce-based sample from Ref. 29. The 
broad lines are visual guides. 


most of the samples, with an increase of the resonance 
frequency when increasing the temperature. The effect is 
largest for the undoped samples, and appears reduced or even 
canceled for single-peaked spectra at high doping, such as 
for the xnqr/xfswt=^-9Q La-based sample. When the 
spectrum features two peaks, the temperature dependence 
of the high-frequency peak seems to be lower than that of 
the low-frequency peak, as can be seen for instance for the 
xnqr/xfswt=0.65 Sm-based sample. An increase of fq 
with increasing temperature is also observed in other iron 
pnictides [45, 46], and goes counter to the expected effect of 
lattice expansion and lattice vibrations [47]. This suggests that 
other changes take place with temperature, for instance intra¬ 
cell atomic displacements, which are beyond the scope of the 
present study and are left as an open question. 

In the following, we focus on the T« 160-300 K range, 
i.e., a temperature range above all low-temperature transi¬ 
tions. Using the data plotted in Fig. 2, the ratio Afq/AT is 
extracted for all spectral peaks of all samples over this temper¬ 
ature range. It is then plotted on Fig. 4 versus the resonance 
frequency of each corresponding peak at room temperature, 
i.e., as far away as possible from low-temperature physics. 
Such an abscissa allows to reflect the frequency branch struc¬ 
ture of Fig. 2 (see dotted boxes on Fig. 4). While the uncer¬ 
tainty on Avq/AT for a given peak (a single point) is large, 
for a given rare earth the peaks belonging to the same branch 
feature a dispersion in Avq/AT which is small enough to 
differentiate the various frequency branches, especially for 
Ce- and Sm-based samples. The tendency is for Avq/AT 


to decrease when going from the L frequency branch to the 
Hi branch (and then to the H 2 /H 3 branches for Sm-based 
samples), whereas for a given branch large differences in the 
doping and low-temperature properties of the samples appear 
to be of lesser importance. This shows that the temperature 
dependence of the electronic charge environment at each ar¬ 
senic site is a local electronic property, rather than an emer¬ 
gent property of the whole sample at a given overall dop¬ 
ing. Note how this is different to previously-obtained 1/Ti 
spin-lattice relaxation rate measurements on La-based sam¬ 
ples [ 21 ], which showed a large doping dependence for both 
L- and H-branch sites. We ascribe this difference to vq being 
sensitive to the charge degree of freedom (a defined configu¬ 
ration of id orbital occupancies) and to 1/Ti probing the spin 
degree of freedom (doping-dependent magnetic fluctuations). 

The fact that the local electronic properties of the undoped, 
magnetically-ordered compound are well-retained across the 
L frequency branch suggests that the low-frequency regions 
may play the key role in establishing orthorhombicity and 
static magnetism, whereas the high-frequency regions may 
play no role at all. To test this hypothesis, the order parameter 
e = [a — h)/{a + b) (a, b: in-plane lattice parameters) of the 
structural transition is plotted as a function of Xnqr/xfswt 
in Fig. 5. In addition to data from Ref. 29 which corresponds 
to our samples, for which xnqr is known, we include several 
other studies from the literature [48-53]. Since the latter do 
not include NQR data, we use the reported phase transition 
temperatures to extract x nqr from our reconstructed phase 
diagram, as shown in Appendix A. The agreement between 
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Figure 4: (Color online) Temperature sensitivity of the quadrupolar frequency, as a function of the peak position at room tempera¬ 
ture. AvqIAT is calculated from measurements at room temperature and at r=150-160 K. The legend indicates xnqr/xfswt ■ For 
xnqr<xfswt, each dotted box groups the points belonging to the same frequency branch (L, Fli, F12, Fla) as defined on Fig. 2. 
For La-based samples, the hatched rectangle about uqKiW.l MFlz indicates the frequency range over which Afq/AT goes from weak 
{xnqr/xfswt = ^-^9) to zero or negative {xnqr/xfswt = ^-9Q). Error bars for the Sm-based samples are not shown, as they exceed the 
actual point dispersion for each frequency branch, suggesting over-estimation due to more complex spectra. 


all studies in Fig. 5 is good, which validates our determina¬ 
tion of the doping. The decrease of e with xmqr is com¬ 
patible with a linear decrease until e=0 for xnqr=xfswt, 
i.e., when the low-frequency spectral weight vanishes. This is 
clear evidence of the low-doping-like regions fully account¬ 
ing for the orthorhombicity of the material, and thus for the 
closely-associated static magnetism. 



Figure 5: (Color online) Order parameter of the structural phase tran¬ 
sition £ = (a — b)/{a + b) (a, b: in-plane lattice parameters) as 
a function of doping. The structural data for all samples is taken 
from Refs. 29, 48-53. The dotted line is a visual guide. Verti¬ 
cal dashed lines indicate the percolation thresholds for 2D and 3D 
nearest-neighbor square-lattice site percolation. 


B. Doping-dependence 

From Fig. 5 and Fig. 8, where all the samples con¬ 
sidered in this article are present, it can be seen that 
xnqr/x_fswt=0.40(4) (Oka-0.03, see Tab. I) is the high¬ 
est doping for which a structural or magnetic transition is 
reported in a La-based sample, and x^qr/xfswt=^-^'^{^) 
(one of our samples) the lowest one for which no such 
transition is observed. This defines a narrow threshold 
xnqr/xfswt~^-T>9{2>) beyond which only superconductiv¬ 
ity occurs. For Ce samples this threshold is about 0.66(6) 
(Zhao-0.1, for which vanishing orthorhombicity is reported). 
For Sm samples it is at 0.67(2), i.e., between 0.75(10) 
(Martinelli-0.2) and 0.65(4) (one of our samples). Such val¬ 
ues agree well with the known percolation thresholds for 
the square-lattice site problem in two and three dimensions, 
which given as dilution thresholds (i.e., the fraction of sites 
to remove) are respectively close to 0.41 and 0.69 [54-56]. 
This strongly suggests that the disappearance of orthorhom¬ 
bicity and static magnetism on doping is ultimately limited by 
the dilution of low-doping-like iron atoms, without influence 
of the superconductivity. Note that the interpretation of the 
above values could be complicated by the fact that each As 
atom probes a square of four Fe atoms, including a minor¬ 
ity of As atoms at the boundaries between LD-like and HD- 
like regions a few nanometers in size (see Sec. V). However, 
the linear variation of Xnqr/xfswt from 0 to 1 upon dop¬ 
ing indicates that it is indeed equal to the relative weight of 
“xNQR=XFSWT”4ike iron atoms. 

The change from two-dimensionality (La) to three- 
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dimensionality (Ce, Sm) indicates an enhanced coupling be¬ 
tween the iron planes when the rare earth is Ce or Sm. Struc¬ 
tural effects can be ruled out since there are only marginal 
differences, especially between La- and Ce-based compounds 
[57]. A more likely explanation is the magnetism of the rare 
earth (Ce, Sm) itself, whose 4/ magnetic moment should in¬ 
crease interplane coupling compared to the nonmagnetic rare 
earth (La). Indeed, magnetic rare earths are known to cou¬ 
ple to the 3d electrons [43, 58-61], and are sufficiently cou¬ 
pled to each other to develop an ordered magnetic state at a 
few kelvins. Note also that they should not significantly alter 
the in-plane couplings, as they lie directly above/below iron 
atoms. Such a situation is compatible with the parent com¬ 
pound having similar Tjy for all rare earths. Indeed, for the 
small ratio r = Jd Jab of the interplane to the intraplane cou¬ 
pling in 1111 compounds [62, 63], Tm is predicted to vary lit¬ 
tle with r. As an example, assuming localized spins 5'= 1/2, a 
tripling of r from 0.001 to 0.003 yields only an increase of Tm 
by 13% [64]. In this picture of increased interplane coupling, 
it is interesting to note that the short-range-ordered (SRO) 
magnetism transition line of Ce and Sm samples shows a 
change in curvature between the dopings 0.34(2) and 0.50(2), 
i.e., about the 2D percolation threshold (see Fig. 3). On cross¬ 
ing the threshold, the in-plane density of low-doping-like re¬ 
gions becomes low enough that in-plane couplings between 
these regions is weakened, effectively increasing the dimen¬ 
sionality. 

In light of the structure, it is expected that in-plane mag¬ 
netic couplings occur between nearest-neighbor (NN) and 
next-nearest-neighbor (NNN) iron sites, i.e., that the J 1 -J 2 
model applies. Square lattice NNh-NNN percolation has been 
argued to occur for La/Sm/Pr-1111 with ruthenium substitu¬ 
tions at iron sites [65-67]. There, spin dilution would cause 
static magnetism to vanish for xiiu^Q.6, close to the theoreti¬ 
cal value 0.593 [68]. However, there is signihcant evidence of 
a nonrandom distribution of Ru and microscopic phase sepa¬ 
ration in Ru-poor and Ru-rich regions [66, 69-73], suggest¬ 
ing that the global Ru content is not the relevant parameter 
to describe percolation. In Appendix B, we propose a simple 
model to reinterpret the NQR data of Sanna et al. [66] by as¬ 
suming a phase separation in Fe-rich and Ru-rich regions. For 
the a;^°'"““*=56% composition, close to which static mag¬ 
netism and superconductivity vanish, we hnd a volume frac¬ 
tion «48% of Ru-rich regions with 80% local Ru concentra¬ 
tion, which is likely too high to allow for any static magnetism 
of the iron atoms embedded in the nonmagnetic Ru matrix 
[74]. Static iron magnetism should instead originate from 
the «52% volume fraction of Fe-rich regions, which harbor 
low-doping-like (LD-like) areas as dehned in the present arti¬ 
cle. Due to the phase separation, the Ru-rich regions should 
be inefficient at preventing percolation of these LD-like ar¬ 
eas. It is thus striking that the fraction of LD-like areas in 
the Fe-rich regions is equal to 0.274(66), i.e., consistent with 
3D NN percolation (threshold equal to 0.31) rather than with 
2D NNh-NNN percolation (0.407). Although further study of 
the role of Ru substitutions is needed, this suggests agreement 
with the behavior observed by us in F-doped samples with a 
magnetic rare earth. Note that comparison with LaFeRuAsO 


using the NQR data of Ref. 75 could not be performed, since 
the very broad magnetic transitions for La-based samples [65- 
67] imply a more complicated distribution of ruthenium in the 
iron planes, for instance a variability of the ruthenium con¬ 
tent in the iron-rich regions of each sample. One possibility 
for the observation of NN percolation is that the elementary 
regions which percolate to yield static magnetism are larger 
than one lattice cell and interact through some effective NN 
interaction, whether local in character or involving the con¬ 
duction electrons. Besides, even in the site-to-site J1-J2 pic¬ 
ture, complications are expected, such as the NNN hoppings 
eventually dominating the NN hoppings at high doping [76]. 


C. Mechanism 

An immediate consequence of the observed percolation be¬ 
havior is that the loss of static magnetism when doping is not 
related to a degradation of the nesting properties of the Fermi 
surface. An even more fundamental question about iron pnic- 
tides is whether the tightly-related structural and magnetic 
transitions are driven by orbital or spin fluctuations [13]. In 
the presence of a magnetic rare earth, not only is 3D mag¬ 
netism obtained on doping beyond the 2D percolation thresh¬ 
old, but orthorhombicity also occurs until the 3D percolation 
threshold (seemingly with an increased dimensionality [29]). 
Considering the crystal structure where R-{0,F)-R (R: rare 
earth) layers separate the FeAs layers, it is hard to explain how 
the unpaired 4/ electrons of the rare earth would promote or¬ 
bital correlations, rather than spin correlations. This suggests 
that the electronic anisotropy of the 1111 system at low dop¬ 
ing is primarily driven by spin fluctuations, eventually yield¬ 
ing the enhanced, slow spin-nematic fluctuations seen below 
the structural transition [77, 78], or at least that there is a de¬ 
cisive feedback of spin fluctuations on orbital fluctuations. 


V. SUPERCONDUCTIVITY 

The high-frequency, high-doping-like regions are natural 
candidates to harbor superconductivity. However, the super¬ 
conducting volume fraction of 1111 materials appears to be 
large even at low doping [26]. Previous spin-lattice relaxation 
measurements also showed that the low- and high-frequency 
regions show a similar decrease of spin excitations below 
Tc [21, 33]. These observations could suggest that the low- 
frequency, low-doping-like regions support superconductivity 
when static magnetism is repressed, which is a crucial point 
with respect to ground state competition and coexistence. 

In La-1111, there is a significant decrease of the superfluid 
density ps below optimal doping [22]. In the case where 
both ground states would coexist, competition for carriers 
would be expected. Such is the case in the 122 compound 
Ba(Fei_a;Coa;) 2 As 2 , where a large enhancement of the pene¬ 
tration depth Xab is observed when the doping is low enough 
for coexistence [79]. is inversely proportional to ps, 
whose decrease is argued to be due to the SDW partially gap¬ 
ping the Fermi surface, i.e., removing carriers. Even though 
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no static magnetism is present in superconducting La-1111, 
it cannot be ruled out that the remaining slow spin fluctua¬ 
tions seen by NMR below optimal doping [30] are related 
to a partial gapping of the Fermi surface over some volume 
fraction of the samples. However, the observed decrease of 
Ps is much faster than that of Tc, yielding a Uemura plot 
(closed orange symbols on Fig. 6 ) where the pnictides show 
no well-defined slope as opposed to underdoped cuprates [80]. 
Taken at face value, this would suggest that superconductivity 
presents large variations of the coupling strength versus dop¬ 
ing. Writing ps=ns/m* with rig the superfluid carrier density 
and m* the effective mass of the Cooper pairs, and consider¬ 
ing the complex multiband structure, an explanation would be 
that changes in m* occur with doping [ 20 ]. 



Figure 6: (Color online) Uemura plot showing the critical tempera¬ 
ture versus All measured data is taken from Refs. 22, 29, 81 (in¬ 
cluding one sample whose Xnqr/xfswt is determined from its 
value, see Appendix A), xnqr/xfswt is shown as a label. The un¬ 
certainty on is given for all samples where uncertainties on Xab 
and xnqr are both available. The dashed line is a visual guide. Dot¬ 
ted lines indicate the behavior observed for hole-doped and electron- 
doped cuprates below or at optimal doping [80, 82]. “s-M” refers 
to the presence of strong static magnetism, and “SC1-I-SC2” to the 
crossover between two superconducting states. 


Here, we argue that the measured superfluid density (usu¬ 
ally using pSR) is reduced for xnqr/xfswt<^ due to the 
proximity effect. Indeed, if only the HD-like regions are in¬ 
trinsically superconducting, then their nanoscale coexistence 
with the LD-like regions implies that the proximity effect may 
play over a large volume fraction, up to the entirety of the LD- 
like regions. Assuming the extremal case, the measured super¬ 
fluid density then corresponds to a 100% SC volume fraction 
originating only from the HD-like regions, i.e.; 


m ^ i Xnqr 
Ps Ps ’ 

XPSWT 


(4) 


where are the measured and intrinsic superfluid den¬ 
sities. Applying this correction to our samples for which 
penetration depth data is available [22, 29, 81], this yields 
a well-defined Tc=f{ns/m*) slope, in-between that of 
hole- and electron-doped cuprates, for the La samples with 


QA9<xnqr/xfswt<^-53 and for the Sm sample with 
xnqr/xfswt=0-9Q (open turquoise symbols and visual 
guide on Fig. 6 ). This single slope indicates a common, well- 
defined superconducting state behavior, suggesting the valid¬ 
ity of this scenario where only the HD-like regions intrinsi¬ 
cally support superconductivity. Even in the absence of static 
magnetism, the LD-like regions do not appear to contribute to 
the condensate. This suggests that in the undoped limit there is 
no underlying superconducting ground state which static mag¬ 
netism would be suppressing. The implications for the mech¬ 
anism responsible for Cooper pairing are also strong, since 
the linearity of Tc on ps would seem to rule out theoretical 
approaches in the weak-coupling limit [83]. 

Note that, for dopings where the material is electroni¬ 
cally inhomogeneous {xnqr/xfswt<^), it cannot be said 
that the observed behavior is that of a material that is sim¬ 
ply less doped than xfswt- The observed Uemura rela¬ 
tion then rather reflects the properties of xpswT-^^^^ su¬ 
perconductivity which is weakened due to the interspersed 
low-doping-like and high-doping-like regions. This seems 
however to be equivalent to a doping variation, since 
the electronically-homogeneous, roughly optimally-doped 
xnqr/xfswt=^-^9 La-based sample obeys the same be¬ 
havior. For La samples with even higher Xnqr/xfswt, 
whose Tc would suggest they are overdoped, the Uemura re¬ 
lation appears to be maintained through the concomitant de¬ 
crease of Tc and ps up to at least xnqr/xfswt=^-33, and 
breaks for xnqr/xfswt=^-9Q which shows a reduced Tc 
and a still moderate superfluid density. This is unlike the be¬ 
havior observed in cuprates, where overdoping tends to re¬ 
sult in a saturation then a decrease of Tc as ps keeps increas¬ 
ing [80]. Here, the Xnqr/xfswt=^-33 La-based sample is 
behaving like an underdoped sample, suggesting that the ob¬ 
served Tc dome is misleading. Note that this is in agreement 
with the study of hydrogen-doped 1111 pnictides by limura 
et al. [44], who report that the decrease of Tc beyond optimal 
doping may be due to a degradation of the nesting of the Fermi 
surface and of the associated spin fluctuations. Concerning 
the xnqr/xfswt=^-90 La-based sample, it is interesting 
to note that the breakdown of the Uemura relation occurs in 
the region of the phase diagram where Tc is depressed (see 
Fig. 3), which corresponds to a crossover between two dif¬ 
ferent superconducting states (SCI and SC2) [44]. Assuming 
that SCI and SC2 may be competing but do not support each 
other, we define in the following 741 ) 1 ( 2 ) j-gspective 

critical temperatures at a given doping, and as their 

respective contributions to the total superfluid density For 
xnqr/xfswt=^-9Q, Fig. 3 indicates Tc^^>Ti^\ While the 
available data is insufficient to conclude, a possible scenario 
is then that the Uemura relation is still obeyed for SCI with 
Tc=Tc^\ and with accounting for the apparent break¬ 
down. For highly-doped Ce-based samples, the breakdown of 
the Uemura relation could have a similar explanation. Finally, 
the insufficient correction of the superfluid density observed 
for the xnqr/xfswt=0-63 Ce-based sample will be linked 
to ground state competition in Sec. VI. 

Our result on the presence of a full-volume proximity effect 
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allows to put an upper boundary on the typical distance 
dn-H separating HD-like regions, with the coherence 
length in the normal material, i.e., the characteristic length 
over which the pair amplitude decays. For similarly poor con¬ 
ductors such as cuprates, can be expected to be on the or¬ 
der of a few nanometers. Since the LD-like regions are unable 
to support superconductivity by themselves, it is possible to 
rule out an enhancement of at low temperature such as that 
observed in Josephson junctions where two superconducting, 
high-Tc cuprate electrodes are separated by an underdoped, 
lower-Tc cuprate barrier in the normal state [84, 85]. The 
small value of dn-H, together with the conducting electronic 
background, favors the strong coupling of the intrinsically- 
superconducting regions [ 86 ]. This may explain why, in the 
absence of static magnetism, a sample with Xnqu/xfswt 
as low as 0.49 (La-based sample) still achieves a Tc as high 
as Ki75% of the optimal value. Such a strong coupling was 
also proposed to explain the resilience of superconductivity 
in Ru-doped Ba-122 [19]. Finally, the high-doping-like re¬ 
gions should be at least of size ^ 5 , the superconducting coher¬ 
ence length. This translates to about 20^0 A in the ab plane 
[87, 88 ]. Therefore, a patchwork of regions a few nanome¬ 
ters across is the most simple arrangement compatible with 
the observed superconducting behavior. A similar conclusion 
was reached by Sanna et al. on the basis of the magnetic be¬ 
havior [4, 32]. 


VI. INTERPLAY OF THE GROUND STATES 

In light of the detrimental effect of iron static magnetism on 
superconductivity (see Sec. HID), ground state coexistence 
appears as a special case where static magnetism from the 
LD-like regions may restrict the spatial extension of super¬ 
conductivity originating from the HD-like regions. Indeed, 
such a behavior seems present in the xnqr/xfswt=Q-^^ 
Ce-based sample. As seen in Fig. 6 , correcting the superfluid 
density of this sample according to Eq. 4 is not enough to 
bring it in agreement with the observed Uemura relation. This 
suggests a loss of superfluid density, on contrary to the sam¬ 
ples without ground-state coexistence. According to /iSR the 
magnetic volume fraction reaches 100 %, with a freezing of 
the iron moments over an extended temperature range («5- 
26 K) [29]. Such a distribution of Tm is a common fea¬ 
ture [4, 6 , 7, 32]. The 100% volume fraction does not mean 
that the whole sample is intrinsically magnetic, which would 
be incompatible with our results, but that the field from the 
magnetically-ordered LD-like regions is also felt in the HD- 
like regions due to the nanoscale separation. At a given Tjv, 
the corresponding magnetic regions and their nanometer-scale 
environment will show up as frozen volume in the pSR mea¬ 
surement. In isovalently-doped Ba(Fei_ 2 ;Rua;) 2 As 2 , the dis¬ 
tribution of Tft was ascribed to a spatial variation of the or¬ 
dered iron moment, with values smaller than about 0.3 ps 
allowing for local coexistence of superconductivity and static 
magnetism [19]. In F-doped 1111 pnictides, the variation of 
the strength of static magnetism could be due to finite-size 
effects and to the disorder inherent to a transition with per- 


colative character (see Sec. IV). The LD-like regions with the 
strongest magnetism (highest Tm) will locally suppress su¬ 
perconductivity by proximity, and possibly also the intrinsic 
superconductivity in the nearest HD-like regions, whereas the 
LD-like regions with weak or zero static magnetism will toler¬ 
ate it. In such a picture some superfluid density is irremedia¬ 
bly lost in or close to the more magnetic regions, in agreement 
with the observed breakdown of the Uemura relation. This 
would also explain the report of a reduced superconducting 
volume fraction in the presence of ground state coexistence 

[4]. 

All our results which pertain to the existence and the in¬ 
terplay of static magnetism and superconductivity are sum¬ 
marized in Fig. 7. In the parent compound, the paramag¬ 
netic state is electronically homogeneous, yielding the well- 
known long-range-ordered magnetic phase at low tempera¬ 
ture. On increasing the doping nonmagnetic HD-like regions 
start to appear, while Tm decreases. The static magnetism 
from the LD-like regions is nonetheless still strong enough 
to suppress any superconductivity that may have emerged 
from the HD-like regions, thus a zero superconducting vol¬ 
ume fraction {Vsc)- For the nonmagnetic rare earth (La), or- 
thorhombicity and static magnetism cannot occur beyond the 
2D nearest-neighbor percolation threshold. The development 
of superconductivity is then no longer hampered anywhere in 
the material, and Vsc jumps to 100 % due to superconduc¬ 
tivity by proximity in the LD-like regions. A homogeneous, 
purely intrinsic superconducting state is eventually reached 
when the LD-like regions have disappeared. For the mag¬ 
netic rare earths (Ce, Sm) static magnetism survives to higher 
doping, allowing for a distribution of Tm that crosses progres¬ 
sively into a low-temperature range; the LD-like regions with 
weaker magnetism are no longer able to suppress supercon¬ 
ductivity originating from the nearby HD-like regions, and 
Vsc starts to grow. This occurs at xmqr/xfswt^O-55, 
i.e., in-between 0.50 (purely magnetic Ce sample) and 0.59 
(magnetic and superconducting Sm sample), with the approx¬ 
imation that the boundary is the same for Ce- and Sm-based 
compounds. On approaching the 3D nearest-neighbor perco¬ 
lation threshold {xmqrJ XFSWT~f>-&^), Tm should be sys¬ 
tematically low before actually vanishing. Intrinsic super¬ 
conductivity in the HD-like regions is then never suppressed, 
nor is superconductivity by proximity in the LD-like regions. 
Again, a homogeneous superconducting state is reached for 
xnqr=xfswt- While this scenario accounts well for the 
qualitative difference between the phase diagrams of non¬ 
magnetic and magnetic rare earths [3, 4, 22, 32], note that 
a small coexistence region cannot be ruled out for La-based 
samples, as suggested by measurements on hydrogen-doped 
samples [89]. Any such coexistence region will be limited to 
the very narrow doping range (xmqr/xfswt=0.39±0.03, 
see Sec. IV) over which static magnetism could be weak¬ 
ened enough to allow for the presence of superconductivity. 
It should also be noted that for Ce- and Sm-based samples at 
intermediate doping, the detailed evolution of Vsc should de¬ 
pend not only on the distribution in the strength of the static 
magnetism but also on the microscopic arrangement of LD- 
like and HD-like regions, which will affect how much super- 
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Figure 7: (Color online) Summary of the interplay of ground states in the doping range where LD-like (hlue) and HD-like (red) regions 
coexist at the nanoscale. The LD-like regions can feature strong (“s-M”) or weak (“w-M”) static magnetism, as well as superconductivity 
through proximity (“p-SC”). The HD-like regions feature intrinsic superconductivity (“i-SC”), which can he suppressed (“i-SC: lost”). Vsc 
is the superconducting volume fraction, while “LRO” and “PT” stand for “long-range ordered” and “percolation threshold”. For the magnetic 
rare earth case, the intermediate doping region (0<l/sc<100%) is simplified by not making a distinction between strong static magnetism 
preventing superconductivity by proximity and destroying nearby intrinsic superconductivity, and by not detailing the behavior close to the 3D 
percolation threshold, where all of the static magnetism is expected to be weak enough for Vsc~i00%. 


conductivity is suppressed. 

VII. CONCLUSION 

The phase diagram of 1111 i?FeAsOi_a;F 2 ; (ii=La, Ce, 
Sm) iron pnictides has been constructed as a function of the 
NQR spectral properties in the paramagnetic state, i.e., the 
properties of two types of local charge distributions associ¬ 
ated to low-doping-like (LD-like) and high-doping-like (HD- 
like) regions. The combination of a local probe technique and 
relative intensity measurements results in high accuracy and 
consistency across all measured samples. This allows to show 
without any ambiguity that magnetic rare earths (Ce, Sm) pro¬ 
mote iron static magnetism to higher doping, with a detrimen¬ 
tal effect on superconductivity. Besides, our approach allows 
for successful inclusion of NQR data from other groups with 
different sample sources, yielding to our knowledge the most 
extensive phase diagram for fluorine-doped La-based com¬ 
pounds, including the upturn of the superconducting transition 
temperature at high doping. 

Using the NQR-defined doping, as directly measured for 
our samples or as derived from our phase diagram for stud¬ 
ies using other techniques, we have investigated the spatial 
origin of static magnetism and superconductivity. It is found 
that the LD-like regions are closely associated to the devel¬ 
opment of orthorhombicity and static magnetism, with the 
upper doping limit set by dilution effects: 2D (respectively 
3D) nearest-neighbor square lattice site percolation is at play 
when the rare earth is nonmagnetic (respectively magnetic). 
The LD-like regions are not intrinsically supportive of super¬ 
conductivity but can harbor superconductivity by proximity, 
originating from the nearby HD-like regions, whenever static 
magnetism is weak enough. In the end, the interplay of ground 
states in 1111 pnictides appears to be well described as a spa¬ 
tial competition between nanoscopic regions, which compete 
to establish the ground state through suppression of supercon¬ 
ductivity by static magnetism, and extension of superconduc¬ 


tivity by proximity effect. 

Our conclusions are compatible with various observations 
taken from the literature, whether associated to static mag¬ 
netism, such as a pseudogap-like feature that survives up to 
TwITO K in superconducting samples and the Nernst signa¬ 
ture of a possible SDW precursor in an underdoped supercon¬ 
ducting sample without static magnetism [90, 91], or associ¬ 
ated to superconductivity, such as the presence of additional 
pinning due to nanoscale regions with suppressed supercon¬ 
ductivity and the possibly-related enhancement of the upper 
critical field Hc 2 and dHc. 2 /dT{Tc) at low doping [92, 93]. 
The well-defined Uemura relation is also in agreement with 
work supporting the strong coupling nature of superconduc¬ 
tivity in iron pnictides [94]. Besides, our scenario eliminates 
the controversy surrounding the “first-order-like” transition in 
the i?=La phase diagram from magnetically-ordered samples 
to superconducting samples. Rather than a smooth evolution 
of the ground state throughout the whole sample volume, the 
conditions for magnetic order or superconductivity are already 
met on a local scale, with the switching of the ground state be¬ 
ing dependent on percolation. 

Through a better understanding of ground state competition 
and of the nature of superconductivity, our results should help 
bring insight on the origin of unconventional superconductiv¬ 
ity in iron pnictides. Further work should address the origin 
and the details of the nanoscopic arrangement of magnetism- 
prone and superconductivity-prone regions, for which some 
theory leads exist [95, 96], as well as the applicability of our 
picture of ground state interplay to other iron pnictide fami¬ 
lies. 
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Rare earth 

Source - first author 

X 

xnqr/xfswt 

La 

[22] - Luetkens 

0.125 

1.53(6) 


[48] - Huang 

0.03 

0.15(7) 



0.05 

0.36(4) 


[49] - Qureshi 

0.045 

0 .22(6) 

Ce 

[50] - Zhao 

0.02 

0.10(7) 



0.04 

0.28(5) 



0.06 

0.51(4) 



0.08 

0.59(3) 



0.1 

0 .66(6) 

Sm 

[53] - Martinelli 

0.05 

0.30(5) 



0.075 

0.35(5) 



0.1 

0.64(6) 



0.2 

0.75(10) 


[51] - Margadonna 

0.05 

0.40(8) 



0.1 

0.52(4) 



0.12 

0.55(4) 


Table II: Samples from the literature which are included in Fig. 5 and 
Fig. 6. X indicates the fluorine doping as given in the source article, 
which is usually the nominal value, xnqr/xfswt is derived from 
Fig. 8. 
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Appendix A: Doping determinations nsing the reconstructed 
phase diagram 


Using the phase diagram obtained in Sec. Ill, it is possi¬ 
ble to extract xnqr/xfswt for samples from the literature 
whose phase transition temperatures are known, without the 
need for NQR data. As shown on Fig. 8, this has been done 
for a number of samples (listed in Tab. II) whose properties 
(orthorhombicity, penetration depth) are of interest for Sec. IV 
and Sec. V. Agreement could be achieved for samples show¬ 
ing more than one transition, with the exception of the sam¬ 
ples from Ref. 53. The latter feature structural phase transition 
temperatures which are unusually high, due to a specific anal¬ 
ysis of synchrotron data which is sensitive to incipient order. 
For these samples, priority in setting x^qr is then given to 
the magnetic and superconducting phase transitions. 


P 

Experiment 

Binomial 


Model 


Wp 

Wp 


L 

0.142(26) 

0.037 

0.274(66) 

0 .001(1) 

0.143(26) 

1 

0.217(26) 

0.191 

0.419(10) 

0.024(13) 

0.229(20) 

2 

0.251(26) 

0.364 

0.240(47) 

0.149(45) 

0.196(24) 

3 

0.189(26) 

0.309 

0.061(25) 

0.408(22) 

0.228(26) 

4 

0.201(26) 

0.098 

0.006(4) 

0.417(81) 

0.204(26) 


Table III: Spectral weights for the local environments seen by NQR 
in SmFe0.44Ru0.56AsO0.85F0.15 (see text). The experimental data is 
taken from Ref. 66, and the experimental uncertainty is assumed to 
be the point size of the original figure. 


Appendix B: Phase separation in Ru-substituted compounds 


Following the observation of microscopic phase separa¬ 
tion in Fe-rich and Ru-rich regions in 1111 samples where 
iron has been partially substituted with ruthenium [66, 69- 
73], we propose a simple model to reinterpret the data of 
Ref. 66. Applying NQR to SmFei-j^RuajAsOo.ssFo.is, Sanna 
et al. obtain spectra featuring low- and high-frequency com¬ 
ponents like ours but also four more peaks corresponding to 
As nuclei with 1/2/3/4 Ru nearest-neighbors (NN). We write 
the corresponding experimental spectral weights as with 
p=F/H/l/2/3/4. For the a:^°'”™“*=0.56 compound, in which 
static magnetism and superconductivity are on the verge of 
vanishing, NQR yields the weights given in Tab. Ill (as well 
as u>“^=0). As then noted there is a tendency towards Ru 
clustering, with being twice larger than expected for a 
binomial distribution (see iCp™ in Tab. III). 

Beyond phase separation into Fe-rich and Ru-rich regions, 
our model is built on these hypotheses: (i) each type of region 
features binomial statistics corresponding to the local Ru con¬ 
centration (ii) the phase separation is sufficiently marked that 
the Fe-rich regions harbor almost all As with zero or one Ru 
NN, and the Ru-rich regions almost all As with 4 Ru NN. In 
the following, all quantities related to the Fe- and Ru-rich re¬ 
gions are respectively indicated by the exponents “i” and “r”. 

We first derive all compositional quantities. The Ru content 
in the Fe-rich regions x\^ is given by: 

where the two equalities derive from hypotheses (i) and (ii). 
This yields =0.276(44), thus the weights in Tab. III. 
The volume fraction of Fe-rich regions is given by: 




(h) w 


exp 

L 


0.518(69). 


The Ru content in the Ru-rich regions is given by: 


(a; 


.4 (i) r 
) = 


Ru, -4 

yielding =0.804(39), thus the weights in Tab. III. The 
calculated Ru content of the whole sample is given by: 

= V^x),^ + (1 - V^)x^Fu = 0.531(22). 

















14 




150 


OJ 100 

3 

H—< 

cc 

Urn 

(D 

Q. 

E 

(U 

H 

50 


I—O— ^ 

y •• 




\sRo(Ce/Sm) 


N,LRO 




HCH 




hO^ 


^P“ 
K>H 




. (Ce/Sm) " ^ 




HCH 



La 

Ce 

Sm 

■ 

our samples 

■ 

• 

A 

- 

literature 

4> 

€ 

V 

“ 


e* doping 
F content 


(Ce - scaled / Sm) 

' t 

♦ ♦ ♦ 

T (La - scaled) 


■ ‘ ■ ' ■ ■ « ■ I ■ * ‘ ‘ < 


0.0 


0.5 


1.0 


1.5 

V /y 
NQR FSWT 


2.0 


2.5 


Figure 8: (Color online) Phase diagram as used for the determination of xnqr/xfswt for samples from the literature, whose phase transition 
temperatures are known hut for which no NQR data is available (open symbols). These samples are listed in Tab. II. The closed symbols 
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The spectral weights for the whole sample are given by: 

= v^w; + (1 - v^)w;, 

with p=L/l/2/3/4, yielding the weights in Tab. III. 

We then check the consistency of the calculation. The con¬ 
sistency of hypothesis (ii) is shown by the following popula¬ 
tion ratios being small: 


w\y^ 


0.005(5), 


u;[(l - V^) 
wlV^ 


0.054(41), 


wlV^ 

wl{l - V-) 


0.015(11). 


This is also reflected in for p=L/l/4, and rea¬ 

sonable agreement is found forp=2/3. The small differences 
should be attributed to the phase boundaries and to minor de¬ 
viations from randomness inside the Fe-rich and Ru-rich re¬ 
gions. Finally, the slightly-reduced value of compared 
to jg jjj agreement with experimental determinations 

[75, 97]. Therefore, our simple model appears to be in good 
quantitative agreement with the experiment. It is likely com¬ 
plementary to that of Ref. 73, which cannot easily distinguish 
the quantitative properties of the iron-rich and rathenium-rich 
regions, but provides microscopic arrangements compatible 
with NQR measurements. 
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